Astronomy & Astrophysics manuscript no. dications 


February 5, 2008 


(DOI: will be inserted by hand later) 





Theoretical evaluation of PAH dication properties 



G. Malloci 12 , C. Joblin 1 , and G. Mulas 2 



O 
O 

(N 

<D 

in 

(N 



oo 

O 
O 

Oh 
6 



3 



Centre d'Etude Spatiale des Rayonnements, CNRS et Universite Paul Sabatier Toulouse 3, Observatoire Midi- 
Pyrenees, 9 Avenue du Colonel Roche, 31028 Toulouse cedex 04 (France) — e-mail: [giuliano .malloci ; 
Christine .joblin] @cesr . fr 
2 Istituto Nazionale di Astrofisica - Osservatorio Astronomico di Cagliari, Strada n.54, Loc. Poggio dei Pini, 1-09012 Capoterra 
(CA) (Italy) — e-mail: [gmalloci ; gmulas] @ca . astro . it 

Received 18 July 2006; 25 September 2006 

Abstract. Aims. We present a systematic theoretical study on 40 polycyclic aromatic hydrocarbons dications (PAHs ++ ) 
containing up to 66 carbon atoms. 

Methods. We performed our calculations using well established quantum-chemical techniques in the framework of the density 
functional theory (DFT) to obtain the electronic ground-state properties, and of the time-dependent DFT (TD-DFT) to 
evaluate the excited-state properties. 

Results. For each PAH ++ considered we computed the absolute visible-UV photo-absorption cross-section up to about 30 eV. 
We also evaluated their vibrational properties and compared them to those of the corresponding neutral and singly-ionised 
species. We estimated the adiabatic and vertical second ionisation energy A/ through total energy differences. 
Conclusions. The A/ values obtained fall in the energy range 8-13 eV, confirming that PAHs could reach the doubly-ionised 
state in HI regions. The total integrated IR absorption cross-sections show a marked increase upon ionization, being on the 
average about two and five times larger for PAHs ++ than for PAHs + and PAHs, respectively. The visible-UV photo-absorption 
cross-sections for the 0, + 1 and +2 charge-states show comparable features but PAHs ++ are found to absorb slightly less than 
their parent neutrals and singly ionized species between ~ 7 and ~ 12 eV. Combining these pieces of information we found that 
PAHs ++ should actually be more stable against photodissociation than PAHs and PAHs + , if dissociation tresholds are nearly 
unchanged by ionization. 

Key words, astrochemistry — molecular data — molecular processes — ISM: molecules — ultraviolet: ISM — infrared: ISM 



^ ■ 1. Introduction 



Several spectroscopic features, observed in absorption or in 
emission in a large number of interstellar environments, are 
hypothesised to provide strong evi dence for the presence 
of large carbon-based molecules ( Ehrenfreund & Charnlev 
2000). Among them are the diffuse inters tellar bands (DIBs, 
Herbid Il995t iFulara & Krelowskil bOOfj . seen in absorp- 



Isolated gas-phase PAHs absorb efficiently visible-UV radia- 
tion and are ther efore lik ely to contribute to the DIBs spectrum 



tion ana are tneretore likely to contribute to me DIBs spectrum 
(Leeer & d'Hendecourl Il985l Ivan der Zwet & Allamandolal 
[Crawford et alJll985h and to the UV-bump and far-UV 



tion in the visible-UV; the arom atic infrared bands (AIBs, 
Peete rs et al ] l2004l lTielensll2005l) . dominating the near and 
mid-IR spectrum of many i nterstella r sources; the e xtended 
red emission (ERE. IWitt & Viihl2004tlwitt et al lgOOfjl) and the 
blue-luminescence (BL, Viih et al. 2004, 2005a b), wide emis- 
sion bands peaking in the red and blue parts of the visible spec- 
trum, respectively; the UV-bump at ~ 220 nm a nd the far-UV 
rise in the general interstellar extinction curve fc)rainel l2003: 
lFitzpatrickl2004l) . 

Based on their spectral properties, their high photostabil- 
ity and their organic nature, polycyclic aromatic hydrocar- 
bons (PAHs) are believed to be abundant in the interstella r 
medium (ISM. IPueet & LegeJl989tlAllamandola et alJl9 89). 



1985; 

rise of the extinction curve d.ToblirJll992t I.Toblin et alJll99l . 
After photon absorption, they undergo a fast radiation- 
less transition (internal conversion or inter-system crossing), 
and subsequently relax radiatively by electronic fluorescence 
and/or phosphorescence (possibly contributing to BL and 
ERE) and IR emission in their active C-C and C-H vibra- 
tional modes (hence pr oducing AIBs, iLeger &Pugell fl984: 
lAllamandola et alJl985h . 

Interstellar PAHs are expected to exist in different charge 
and hydrogenation states depending on the physical conditions 
of the host environment llTielensl2 005). In particular, the possi- 
ble prese nce of PAH di cations (PAHs ++ ) in the ISM, first pro- 
posed by Leach ( 1986), has recently received further interest, 
follo wing the proposa l that PAHs ++ could be the carriers of 
ERE j Witt et al.l2006l). Theoretical works on the IR properties 
of a few PAHs ++ (Elli nger et al.lll999t iBauschlicher & Bakes 
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2000: IPauzat & Ellin getfcool show that their vibrational pat- 
terns are similar to those of their parent singly ionized species 
(PAHs + ). From modelling studies it is predicted that spe- 
cific values of electron density and radiation field intensity 
ought t o lead to a population of PAHs ++ in excess of that of 
PAHs + jBakes et alJkoOlalbl) . Therefore PAHs ++ would also 
contribute to the AIBs. 

While the first ionisation energies 7 + for more than 
200 PA Hs containing up to 48 c arbon atoms are well docu- 
mented (Eilfeld & Schmidt 1981), only a few experimental de- 
terminations of the double ionisation energy 7 ++ of some PAHs 
have been published to date. Using photon impact, el ectron 
impact and charge stripping techniques, Tobi ta et all (^994) 
performed a comprehensive study of 7 + and 7 ++ for 21 dif- 
ferent PAHs up to the size of coronene (€24^2). More re- 
cently, Schrod er et alJ J200ll) presented a combination between 
photo-ionisation, charge stripping, reactivity studies and quan- 
tum-chemical calculations to obtain the single and double ion- 
isation energies of corannulene (C20H10) and coronene. 

From the known values of 7 + and 7 ++ , it was suggested 
jLeachll 19861 Il987l Il989l Il995l Il996h that PAHs ++ could be 
formed in diffuse HI regions. The proposed formation mecha- 
nism is a sequential two-photons absorption process in which 
the ionisation of a PAH is followed by a further ionisation of 
PAH + . The minimum energy needed for the first step is the sin- 
gle ionisation energy 7 + , while the second step requires a pho- 
ton energy of at least A7 = 7 ++ - 7 + ; the maximum amount 
of internal energy in the newly formed PAHs ++ is therefore 
13.6 eV - A7. Since the total energy of the dication is larger 
than the sum of the total energie s of two singly io nised frag- 
ments, PAHs ++ are metastable jRosi et aljEiol and could 

be dissociated via Coulomb explosion as shown by exper- 

1 1 I I r J 1 

each et alJ I1989albl) . However, theoretical calcula- 
tions on th e fragmentation pathways of the benzene dication 
jRosi et alJl2004l) show that this charge separation mechanism 
should not be a significant channel in most astrophysical en- 
vironments. This is consistent with the observed production of 
significant amounts of PAHs ++ in soft ionisation experiments 
jLedinghametalJfl 999). Therefore PAHs ++ could be found in 
the ISM. 

To evaluate t he possible contribution of PAHs ++ to ERE 
IWitt et all J2006I) computed the first low-lying electronic tran- 
sitions of five representative PAHs ++ up to the size of ovalene 
(C32H14), showing that they fall indeed in the spectral range 
of the observed ERE. Laboratory measurements of the yield of 
optical fluorescence by PAHs ++ are however needed to test this 

hypothesis. To the best of our knowledge, with the exception 

I 1 1 "i 

of the work by Witt et al. (2006), no other electronic spectra of 

PAHs ++ have been published to date. In particular, a detailed 
study of the electronic excitation properties for a large sam- 
ple of PAHs ++ up to the energies excitable in typical interstel- 
lar sources, together with their vibrational analyses, has been 
missing until now. 

In the last few years we have started a long-term project 
to produce an atlas of synthetic absorption and emission spec- 
tra of spec ific PAHs, to be compared with astronomical ob- 
serva t ions Jjoblin et alJl2002l: iMalloci et ai1l2003l: iMulas et all 
l2003H2006bycir This requires the use of some key molecu- 



lar parameters of PAHs as a basis to run Monte-Ca rlo simu- 
lations jBarkerll 19831 lMulaslfl998l Ijoblin et ai]l2002l) of their 
photophysics in the ISM. We report here these fundamental 
data for a sample of 40 PAHs ++ ranging in size from naphtha- 
lene and azulene (CioHs) to circumovalene (C66H2o). Even if 
PAHs and related species containing from less than 50 to more 
than 200 carbon atoms are expected to be present in the ISM 
(Boulanger 1999), we here restricted ourselves to species con- 
taining up to a maximum of 66 carbon atoms, since computa- 
tional costs scale steeply with dimensions. Follo wing our previ- 
ous work on neutral, cationic and anionic PAHs ( Ma lloci et alJ 
120041 l2005t) . we computed the absolute photo-absorption 
cross-section for the doubly-ionised species. For these same 
molecules we also evaluated the IR spectral properties extend- 
ing the sample of specific PAH^ + for which they were pre- 
vious l y available (|EJhngere^alj fl99^: IBauschlicher & B akes 
2000| IPauzat & Ellingeg|2002K New calculations, together 
with our previous results at the same level of theory for the 
corresponding PAHs and PAHs + , enabled us to estimate, via 
total energy differences, the adiabatic and vertical 7 + and 7 ++ . 
Although the computational approach we used for this purpose 
is not sufficient to predict absolute ionisation energies within 
chemical accuracy (+0.1 eV), it is in reasonable agreement 
with the available experimental data and can be easily extended 
to larger PAHs in the near future. 

Section |2] contains the technical details for the calcula- 
tion of the ionisation energies and the vibrational properties 
(Sect. l2.lt and the photo-absorption spectra (Sect. 12.2V The re- 
sults obtained are presented in Sect.[3]and discussed in Sect.|4] 
Our conclusions are reported in Sect.|3] 

2. Computational details 

Due to the large number of electrons in the molecules con- 
sidered (416 in the largest one, circumovalene), an ab ini- 
tio study based on the direct solution of the many-electron 
Schrodinger equation currently has prohibitive computational 
costs. We here used the density functional theory (DFT, 
Jones & Gunnarsson 1989) and its time-dependent extension 
(TD-DFT, Maraues & Gross 2004), which are the methods of 
choice for the study, respectively, of the ground-state and the 
excited-state properties of such complex molecules as PAHs. 

2.1. Ionisation energies and vibrational properties 

Adiabatic single and double ionisation energies (AIEs) are 
evaluated as the difference between the total energies of the 
neutral and the corresponding cation and dication in their most 
stable configurations, respectively. These quantities, therefore, 
take into account the structural relaxation of the molecule after 
each ionisation process. The single and double vertical ionisa- 
tion energies (VIEs), computed at the optimised geometry of 
the neutral molecule, neglect structural relaxation but include 
wave-functions relaxation following the removal of one and 
two electrons. To estimate AIEs and VIEs we used the hybrid 
B3LYP functional, a combination of Becke's three parameters 
exchange functional jBeckel fl993) and of the Lee- Yang-Parr 
gradient-corrected correlation functional ( Lee et alJ 1988). 
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Hybrid DFT functionals are a little more expensive than other 
exchange-correlation functional but yield better results. The 
B3LYP functional, in partic ular, is widely used in the study of 
PAHs and related species (Martin et al. 1996; Langhoff 199' 
Wiberdl997tlBauschlicher & Langhoffll997tlKato et all 199' 



Kato & Yamabe 2002; Dessent 2000; Bauschlicher & Bakes 
20001 l200lt ISchroder et alJ 1200 it iRienstra-Kiracofe et al 






Rosi et al. 


2004; 


Woon & Park 


Witt eiaU t 


>006; 


Kadantsev et al. 



To assess the impact of this choice we also consid- 
ered some more recently de veloped functio nals, such as the 
so-ca lled B-97 JBeckdll997l) and HCTH jHamprecht et alJ 



1998), but no significants improvements, if any, were obtained 
with respect to B3LYP for the calculation of the ionisation en- 
ergies of PAHs. 

For this part of the work we used the Gaussian -based DFT 
module of the NWChem code ( Ap ra et all 120051) . Geometry 
optimisations were performed using the relatively inexpen- 
sive 4-3 1G basis set, followed by the full vibrational analy- 
ses to confirm the geometries obtained to be global minima on 
the potential energy surface and to evaluate zero-point-energy 
(ZPE) corrections. In the study of the vibrational proper- 
ties of PAHs, the combination B3LYP/4-31G was proven 



scale-factor jLanshofll 


1996; Bauschlicher & Langhoff 


19971 


iBauschlicher & Bakes 


2000, 1200 it 


Bauschlicher 2002 


). We 



thus adopted the same scaling procedure derived by these au- 
thors to obtain the IR spectra of all the PAHs under study in 
their 0, +1 and +2 charge states. Since for open-shell systems 
analytical second derivatives are not yet implemented in the 
currently available version of NWChem (4.7), we performed 
the vibrational analyses for all the cations using the Gaussian03 
quantum chemistry package (Frisch et al. 2003). 

We then started from the optimal geometries and the cor- 
responding self-consistent-field solution obtained in the pre- 
vious step to refine the optimisation with the 6-3 1+G* ba- 
sis, a valence double zeta set augmented with d polarisation 
functions and s a nd p diffuse functions for each carbon atom 
jFrish et alJll984l) . Although basis set convergence is not yet 
expected at this level, in view of the large systems under study 
some compromise between accuracy and computational costs 
had to be made. Thus, while the B3LYP/6-3 1 +G* level of the- 
ory gives a good a greement with experiment for the e lectron 
affinities of PAHs (Dessent l2000t iMalloci et all 1 2005). from 
previous works {S chroder et al. 2001; Deleu ze et alJ 12003: 
IWoon & Parkll2004l:lKadantsev et alJl2006l) it is clear that this 
same level is certainly unable to predict absolute ionisation en- 
ergies within chemical accuracy (+0.1eV). For example, using 
the 6-31 lg** basis set, Schrode r et alJ J200lh obtained mean 
deviations of -0.3 and -0.8 eV by comparing their computed 
values of the single and double AIEs of corannulene (C20H10) 
and coronene (C24H12) with their measured photoionisation 
thresholds. Analogously, the first AIE of five PAHs in the 
oligoacenes series from naphthale ne to hexacene (C4„ + ?H 2 B +4 
with n=2, 3, 4, 5, 6), computed by Kadan tsev et alJ J20061) us- 
ing an even larger basis set, underestimate the experimental val- 
ues by 0.2-0.4 eV; the authors concluded that a better descrip- 



Table 1. Effect of basis set incompleteness on the evaluation 
of single and double AIEs of naphthalene and fluorene at the 
B3LYP level of theory. The reported values (in eV) do not in- 
clude ZPE corrections. 



Basis 


Naphthalene (C 10 H 8 ) 


Fluorene (C^H^) 


set 


single IE 


double IE 


single IE 


double IE 


4-3 lg 


7.61 


20.75 


7.38 


19.72 


6-3 1+g* 


7.80 fli 


20.99* 


7.56 


19.96 


cc-pvdz 


7.75 1 ' 


20.95 


7.51 


19.89 


cc-pvtz 


7.85'' 


21.09 


7.61 


20.04 


Exp. 


8.144±0.001'' 


2\.5±0.2 e 


7.91±0.02 rf 


21.0±0.2 e 



a The use of the 6-3 1+G** basis set (including polarisation functions 
on both C and H atoms) leads to the same value iWoon & Parkl2 004): 
b Using the 6-31 1G** basis set (which is supplemented with a 
third layer of valence functions and includes polarisation functions) 
Schro der et ailfoOOlh obtain 7.82 and 21 .04, respectively; 
c At the same level, Deleuze et al. 1 2003) report 7.75 and 7.83 e V; 
d Recommended values (NIST Chemistry WebBook, Lias 2005 
c Derived from photon impact experiments ( Tobi ta et aljl994 



tion of electron correlation is needed to reproduce experimental 
IEs. 

To estimate the effect of basis set incompleteness, we per- 
formed some reference computations for the AIEs of naph- 
thalene (CioHg) and fluorene (C^H^), using the correlation 
consistent polari sed valence basis sets cc-pVDZ and cc-pVTZ 
3G989). These bases have been used in a bench- 



unmng 



|d 

mark-quality study of the ionisation energies of benzene and 
oligoacenes, to obta in the impressive accuracy of 0.02-0.07 eV 
(foeleuze et al.l F2003). Even if the cc-pVDZ set is less accu- 
rate than the 6-3 1 +G* set and the cc-pVTZ is far from com- 
plete, this comparison provides us some insight into basis set 
requirements. As to zero-point-energy corrections, these same 
authors found that typical values for the PAHs considered are 
in the range 0.01-0.03. Our B3LYP/4-31G values were com- 
puted to be of the same order of magnitude and were therefore 
omitted in the evaluation of the IEs, since they are smaller than 
the accuracy of the method use d. Although the basis s et limit 
is not reached in Tabled (e.g., Kada ntsev et alJEo06l report 
7.88 eV for the first AIE of naphthalene), it is clearly seen that 
the performance of the 6-3 1 +G* basis is intermediate between 
those of Dunning's two basis sets. We thus expect our theo- 
retical predictions to systematically underestimate single and 
double AIEs roughly by 4 and 5 % respectively, which is suffi- 
ciently accurate for the purposes of this work. 

The optimised ground-state structures were all planar, 
with the exception of acenaphthene, fluorene and corannu- 
lene. Whenever possible, molecular symmetry was assumed 
during the calculations. In the case of the more symmetric 
molecules considered, symmetry breaking was observed upon 
single and double i onisation, as expe cted from Jahn-Teller dis- 
tortion JToriill999l:lkato et al.ll999l) . More specifically, we ob- 
tained a symmetry reduction from Dy, to D2h for coronene, 
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Table 2. 

work. 



Complete set of different PAHs considered in this the photo-absorption cross-section cr(E) of neutral and 

singly-ionised PAHs using the real-time real-s pace TD-DFT 
imple mentation of the octopus computer code dMarques et alJ 
2003). We performed these calculations using the widely 
used local-density approximation (LDA), specifically with 
the exchange-corre lation energy density of the homoge- 
neous electron ga s (Ceperlev & Alder 1980) parametrised by 
IPerdew & Zuneed dl98lh . The comparison wit h the e xperi- 
mental data avail able for a few neutral PAHs J.Tohlinl [i992: 
LToblin et alJ fT992) shows that our results reproduce the over- 
all far-UV behaviour up to about 30 eV. This includes the 
broad absorption peak dominated by <x* <— cr transitions, 
which matches well both in position and width. Concerning 
the low-lying excited states of n* <— n character occurring 
in the near-IR, visible and near-UV spectral range, the com- 
puted vertical excitation energies are precise to within a few 
tenths of eV (in the range 0.1-0.4 eV) which are indeed the 
typical accuracies achievable by TD-DFT using h ybrid or gra- 
dient-correc ted exchange-correlation functionals ( H irata et alJ 
Il999l l2003h . We proved these data to be reliable enough for 
detailed modelling of PAHs in different charge states when lab- 
oratory data are missing (Mulas et al. 2006a b c). Starting from 
the ground-state geometries obtained as described in the previ- 
ous section, we extended here similar calculations to PAHs ++ . 

In the scheme we used the wave-functions are represented 
by their discretised values on a uniform spatial grid. The 
molecule is perturbed by an impulsive electromagnetic field 
and the time-dependent Kohn-Sham equations are solved in 
real time. From the knowledge of the time-dependent induced 
dipole moment, cr(E) directly follows by Fourier transform 1 . 
The c alculations we perf ormed for PAHs ++ are almost the same 
as in iMalloci et alJ d2004l) . where a t horough description can 
be found. As already done for anions (Mall oci et al.l2005l) . we 
added an absorbing boundary, which quenches spurious res- 
onances due to standing waves in the finite simulation box 
dYabana & Bertschll999HMaraues et all2003l) . 

The photo-ab sorption spectra of a few PAHs ++ reported by 
IWitt et al] fcOOdl were obtained with a different implementa- 
tion of TD-DFT, which is based on the id entification of the 
poles of the linear response function (Casida 1995). In this lat- 
ter scheme, computational costs scale steeply with the number 
of required transitions, and the excitation sp ectra are limited to 
the lo w-energy part of the spectrum (e.g., iHirata et"aflll999l 
2003). The a dvantages of the real-time propa gation method 
used here are discussed by [Lopez et"afl d2005l) . From the as- 
trophysical point of view, the main step forward we achieved 
using this approach lies in the spectral range covered, that ex- 
tends up to the energies excitable in a typical interstellar source. 
At the same time, the main drawback is that we do not obtain 
independent information for each excited electronic state, such 
as its symmetry and its description in terms of promotion of 
electrons in a one-electron picture. 



Name (Formula) 


Symmetry group 


Azulene (CioHg) 


D 2h 


Naphthalene (CioHg) 


D 2h 


Acenaphthylene (C^Hg) 




Biphenylene (C^Hg) 


D 2h 


Acenaphthene (C12H10) 


C2V 


Fluorene (C13H10) 


C2.V 


Anthracene (C14H10) 


D 2h 


Phenanthrene (C14H10) 


C2V 


Pyrene (Ci 6 H 10 ) 


D 2h 


Tetracene (Ci 8 Hi 2 ) 


D 2h 


Chrysene (C 18 H 12 ) 




Triphenylene (CigH 12 ) 


D 3h 


Benzo[a]anthracene (C18H12) 


C s 


Corannulene (C20H10) 


c 5 , 


Benzo[a]pyrene (C20H12) 


c s 


Benzo[e]pyrene (C20H12) 


C2V 


Perylene (C20H12) 


D 2h 


Anthanthrene (C22H 12 ) 


C211 


Benzo[g,h,i]perylene (C22H12) 


D 2h 


Pentacene (C22H14) 


D 2h 


Coronene (C 2 4H 12 ) 


D 6h 


Dibenzo[b,def]chrysene (C24H14) 


Cib 


Dibenzo[cd,lm]perylene (C26H14) 


D 2 h 


Hexacene (C26H16) 


D 2h 


Bisanthene (C28H14) 


D 2h 


Benzo[a]coronene {Ciziiw) 


D 2h 


Dibenzo[bc,kl]coronene (C30H14) 


D 2h 


Dibenzo[bc,ef]coronene (C30H14) 


C-2v 


Terrylene (C 3 oHi 6 ) 


D 2h 


Ovalene (C32H14) 


D 2h 


Tetrabezo[bc,ef,kl,no]coronene (C36H 16 ) 


D 2h 


Circumbiphenyl (C38H16) 


D 2h 


Circumanthracene (C40H16) 


D 2h 


Quaterrylene (C40H20) 


D 2h 


Circumpyrene (C42H16) 


D 2h 


Hexabenzocoronene (C42H 18 ) 


D 6h 


Dicoronylene (C48H20) 


D 2h 


Pentarylene (CsoHs) 


D 2h 


Circumcoronene (C 54 Hi8) 


C6h 


Circumovalene (C66H20) 


D 2h 



circumcoronene and hexabenzocoronene, from D3h to D2 V for 
triphenylene and from D5 V to C s for corannulene. The complete 
set of the 40 molecules considered and their symmetry point 
groups are reported in Table[2] 

2.2. The real-time TD-DFT method applied to PAHs ++ 

TD-DFT calculations were shown to be a powerful tool to 
calculate electronic excitation properties f o r neutral PA Hs 



2006 


as well as radical ions ui 


3 to large species ( Hirata et al. 


1999 


200U 


alasinski et alJ 


2000; Halasinski et al. 2003; 


Weisman et al. 


200 lL 120031 l2005t IWitt et al] l200d). In our 



previous studies ( Mall oci et al] 120041 120051) we calculated 



1 For further details the reader is referred to the octopus web site 
http : / /www . tddf t . org/programs/octopus . 
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3. Results 

All neutral and singly-ionised species were computed as sin- 
glet and doublet, respectively, while for dications we computed 
both their singlet and triplet ground states. This enabled us to 
predict the relative energies of electronic states having differ- 
ent multiplicities. For almost all of the molecules under study 
(34 out of a total of 40), our calculations predict the dication 
overall ground-state to be the singlet. In the other cases, the di- 
cations of acenaphtylene, triphenylene, corannulene, coronene, 
hexabenzocoronene and circumcoronene, the predicted en- 
ergy difference between the singlet and triplet state is in the 
range 0.03-0.36 eV, which cannot unambiguously identify the 
ground state, given the accuracy of the method we used. In 
this respect DFT-based methods are known to be limited in 
the evaluation of accurate energetic ordering of very close 
electr onic states having different multiplicities JSchroder et alJ 
120011) . In addition, the B3LYP functional is thought to be bi- 
ased towards higher spin states (Bauschlicher & Bakes 2000) 
and better levels of wave function favour the singlet state rela- 
tive to the triplet dRosi et all2004l) . 

3.1. Single and double ionisation energies 

For each species we derived single and double AIEs, respec- 
tively as total energy differences between the geometry opti- 
mised cations and dications and the total energy of the geome- 
try optimised neutral. The corresponding VIEs were evaluated 
as the differences between the single point energy of the cation 
and dication at the optimised neutral geometry and the total en- 
ergy of the optimised neutral. The adiabatic second ionisation 
energy AI is simply given by I^-I^ d while its vertical value 
is obtained through the difference between the single point en- 
ergy of the dication at the optimised cation geometry and the 
total energy of the optimised cation. These results are listed 
in the Table IaTTI of the Appendix and sketched in Fig. ^ as a 
function of the number of carbon atoms in the molecule. The 
differences between vertical and adiabatic ionisation energies 
are found to decrease as the size of the molecu le increases. I n 
all cases however, as already found by Schrode r et alJ (12001 ). 
these values are relatively small, being of the order of 0.1 eV 
for first and second IEs and falling in the range 0.1-0.5 eV for 
double IEs. 

3.2. IR spectral properties 

All the computed harmonic vibrational frequencies of the 
PAHs ++ considered, as well for their corresponding anions, 
neutrals and cations, can be found in our on-line database 
of computed spectr al properties of PAH s, which is presently 
under construction jMalloci et alJl2006i) . We report here only 
comparative plots for the IR properties of the PAHs under 
study in their 0, +1 and +2 charge st a tes. In agreement with 
previous works j Ellinger et alJ 1 19991 iBauschlicher & B akes 
2000; Pauzat & Ellingerl bOOaT we found that while the har- 
monic vibrational frequencies are only slightly affected by sin- 
gle and double ionisation, absolute and relative IR absorp- 
tion cross-sections undergo significant variations with charge. 




6 r 1 

10 20 30 40 50 60 70 
Number of carbon atoms 

Fig. 1. Comparison between experimental (asterisks) and com- 
puted (diamonds) first (a) and second (b) AIEs for the 40 PAHs 
considered in this work as a function of the number of carbon 
atoms i n the molecul e. The continuous line is an empirical es- 
timate llTielenJ2005L p. 191). 



Figure|5]displays the fractions of the total integrated IR absorp- 
tion cross-section in the four spectral ranges 2.5-3.5, 5-10, 
10-15 /mi and > 15 /mi as a function of the total carbon 
atoms in the molecule. These wavelength intervals are some- 
what arbitrary but roughly contain the spectral signatures of, re- 
spectively: C-H stretchings (2.5-3.5 /mi); C-C in-plane stretch- 
ings and C-H in-plane bendings (5-10 /mi); C-H out of plane 
bendings (10-15 /mi); in-plane and out-of-plane cycle defor- 
mations of the whole carbon skeleton (> 15 /mi). The low- 
est panel of Fig. |5] shows the absolute value of the total inte- 
grated IR absorption cross-section divided by the number of 
carbon atoms Nc- In Fig. [2 we omitted the vibrational data 
we obtained for the more symmetric molecules, for which 
Jahn-Teller distortion is observed upon ionisation. For exam- 
ple, in the specific cases of the cations of triphenylene and 
coronene, two different Jahn-Teller distort ed structures are 
known to exist dToriill 19991 llCato et al.lll999i) . These two con- 
figurations are very close in energy and the determination of 
the most stable structure is a delicate problem that may com- 
plicate DFT-ba sed calculations and severely affect the result- 
ing IR spectra JOomens et al.ll200ll) . We report in Fig. [3] the 
corresponding data of Fig.|2]given only for a specific series of 
PAHs, namely the oligorylenes perylene, terrylene, quaterry- 
lene and pentarylene (Cio«H4„ + 4, with n=2, 3, 4, 5). 
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As to the accuracy of our B3LYP/4-31G data, this method 
seems to overest i mate the intensities of C-H stre tching modes 
jLanghofilll99rjl iBauschlicher & Langhofflll997l) . The evalu- 
ation of absolute IR absorption cross-sections is a difficult 
task since they are known to be strongly affected by the 
specific choic e of the basis se t employed llLanghofflll996t 
IBauschlicher & Langhofj|l997l> . This holds true, in particu- 
lar, for open-shell species, that may present problems of spin 
contamination or charge localisation artifacts. As an example, 
the total integrated IR absorption cross-sections we obtain for 
pyrene neutral, cation and dication are 464.4, 784.6 and 1 758.3 
km/mol, while the calculations of Elling eret alJ i 19991) per- 
formed at the higher B3LYP/6-31G* level, give 399.3, 720.9 
and 1613.5 km/mol, respectively. The intensity ratios in the 
same wavelength intervals introduced above are almost coin- 
cident among the two sets of results for Ci6Hiq and Ci6Hiq ++ . 
In the case of the cation Ci6H"io + , instead, these values amount 
to 2, 68, 27 and 3 %, to be compa red with 1, 71, 19 and 9 %, as 
obtained by Ellinge r et alJ Jl999l) . 



3.3. Photo-absorption spectra 



As for the computed vibrational properties, all the theoretical 
photo-absorption spectra of the PAHs ++ considered and their 
neutral and singly-ionised counterparts can be found in our 
on-line database. For the low-lying n* <— n transitions occur- 
ring in the near-IR, visible and near-UV spectral ranges, the 
application of this approach to PAHs and PAHs* was shown to 
be as accurate as previously published theo retical results, com- 
pared to the available experimental data (Mal loci et alj F2004, 
2005). However, while we established the use of the octopus 
code to calculate reasonably accurate photo-absorption spectra 
of neutral PAHs up to the vacuum-UV spectral range, exper- 
imental spectra in the far-UV are needed for a similar direct 
validation for PAH ions. The computed vertical excitation en- 
ergies are expected to be precise to within a few tenths of an 
eV, which are the typical accuracies of TD-DFT calculations 
employing the currently available exchange-correlatio n func- 
tional dHirata et all 1 19991 l2003t IWeisman et aflkOOll 120031 
120051) . More specifically, depending on the character of the spe- 
cific transition, TD-DFT can reach the same accuracy of more 
sophisticated and expensive wave-function based methods or 
show substantial errors, as in the c ase of the lowest s hort-axis 
polarized states in the oligoacenes (Gri mme & Pa rac 2003). 

The spectra obtained for PAHs ++ can be compared with 
the corresponding ones of PAHs and PAHs + . Fig. |4] shows 
weighted averages of the electronic photo-absorption spec- 
tra of the same sample of PAHs, respectively in the neutral, 
cationic and dicationic state. To make apparent the changes 
occurring for specific molecules, Figs. [5] and |6] compare 
the photo-absorption cross-sections of benzo[a]anthracene 
(C18H12) and circumanthracene (C40H16) as a function of in- 
creasing positive charge state. 
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Fig. 2. Comparison between the fraction of the total integrated 
IR absorption cross-section in the spectral ranges 0-5, 5-10, 
10-15 pm and > 15 /im for the neutrals (diamonds), cations 
(crosses) and dications (asterisks) of the PAHs considered 
as a function of the number of carbon atoms Nc- The low- 
est panel shows the absolute total integrated IR absorption 
cross-sections / (in km/mol), divided by Nc- 



4. Discussion 

Concerning the computed ionisation energies, the values we 
found for the second AIEs confirm that PAHs ++ can be ex- 
pected to be produced by two-step ionisation in the ISM. 

As to the IR properties, the positions of the bands 
are relatively similar for the same PAH in different charge 
states, as already noted in previous studies on a few PAH 
dications (Ellinger et al. 1999; Bauschlicher & Bakes 2000; 
IPauzat & EllingeJl2002l) . We confirm that there are system- 
atic shifts in some band positions as a function of ionisation 
state, which only slightly change the general pattern. We here 
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Fig. 3. Same as Fig. |2] restricted to the PAHs of the 
oligoarylenes series, i.e., perylene, terrylene, quaterrylene and 
pentarylene (Cio„H4„ + 4, with n=2, 3, 4, 5). 



concentrate on the much larger variations in the integrated ab- 
sorption cross-sections. The latter change, for some bands, by 
one order of magnitude, due to drastic changes in charge distri- 
bution following ionisation. Our sample of molecules is quite 
heterogeneous and Fig. 13 may be somewhat confusing due to 
crowding. However, general trends versus molecular size and 
versus charge can still be observed, and are more evident upon 
examination of sequences of molecules of the same class, as 
Fig- E] for the oligorylenes series clearly shows. From the top 
panels of Figs.|2]and[3]we see the well-known collapse of the 
intensitie s of C-H str etching modes for singly-ionised species 
foe Frees et alJll993l Eanghofdll996l). The same is observed 
for doubly-ionised species ( Bauschlicher & Bakes 2000) with 
the exception of small PAHs up to ~ 16 carbon atoms, in which 
the inten sity of these modes is partly restored by dou ble ion- 
isation (Ellin ger et aDll999t IPau zat & Ellinger 120021) . A net 
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Fig. 4. Comparison between the weighted sum of the com- 
puted photo-absorption cross-sections <t(E) (expressed in 
Megabarns/Nc, 1 Mb = 10" 18 cm 2 ) for PAH neutrals, the corre- 
sponding cations and dications. We used as weight the inverse 
of the total number of carbon atoms Nc of each molecule. 



increase in the intensity of the in-plane C-H band and C-C 
stretching vibrations is observed for PAHs + with respect to 
neutral species in the 5-10 //m window. This occurs even to a 
larger extent for PAHs ++ , whose absorption spectra tend to be 
dominated only by the features in this region. The opposite is 
observed in the 10-15 fj.m interval, where the out-of-plane C-H 
vibrations are the second major component for neutral species 
but tend to be less important for PAHs + and PAHs ++ . Bands at 
wavelengths larger than 15 /im contribute for less than 10 % of 
the total integrated IR absorption cross-section /. The absolute 
value of the latter quantity, for all of the molecules considered, 
shows a marked increase upon ionisation, as clearly demon- 
strated by the lowest panel of Figs. [2] and [3] On the average we 
found that PAHs ++ absorb about 2.5 + 1 .0 times more than their 
parent singly ionised molecules and about 4.7 + 1.1 times more 
than their parent neutrals. 

With respect to the electronic absorption properties, inspec- 
tion of Figs. |5] and |6] shows that all the spectra for specific 
molecules show comparable features in the UV range, dis- 
playing the^m£_bejiaviour _alreadyknownfor neutral species 
(JoblrnT ll992t Ijoblin et al]ll992t lLeacb1ll995i) . i.e., increas- 
ing smoothly at energies above ~ 7.4 eV to a maximum at 
~ 17.4 eV in a single broad absorption peak composed of 
cr* <— cr, &* <— n, n* <— a and Rydberg spectral transitions. 
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Fig. 5. Photo-absorption cross-section cr(E) of 
benzo[a]anthracene (C18H12) in the three charge states 0, 
+ 1 and +2. Units for cr{E) are Megabarns (1 Mb = 10" 18 cm 2 ). 



The onset of this latter broad absorption moves blueward and 
becomes steeper with increasing positive charge, following the 
increase in subsequent ionisation potentials. This translates in 
a systematic decrease of the cross-section between ~7 and 
~12 eV with increasing positive charge, as shown in Fig. Q 
In the near-IR, visible and near-UV spectral ranges, the dica- 
tion has strong absorption features of n* <— ji character like the 
singly-ionised species. Such a behaviour is an expected conse- 
quence of the different one-electron energy levels and of the 
resulting first few electric dipole-permitted transitions. 

More generally, Fig.@]shows that the three averaged spec- 
tra present very similar features. The main difference is the 
progressive decrease, with increasing positive charge, of the 
cross-section between ~1 and ~12 eV, as expected from the 
behaviour of the individual spectra in this same energy inter- 
val (cf. Fig. 0. Both the former similarity and the latter dif- 
ference have interesting astrophysical implications for the UV 
range of the spectrum: apart from small details (i.e., the pre- 
cise position of single, specific transitions of a single, spe- 
cific molecule), the contribution of any given PAH will de- 
pend only slightly on its ionisation state. The average over a 
sample of molecules has the main effect of smoothing out the 
detailed structure resulting from the contributio n of each sin- 
gle species. This confirms our previous fin dings (Jobli nil 19921 
Ijoblin et alJfl992l iMalloci et al.ll2004 12005I) . extending them 
to PAHs ++ : if any PAHs of small to medium size, be they neu- 



Fig.6. Same as Fig.|5]for circumanthracene (C40H16). 
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Fig. 7. Comparison between the integrated values between ~7 
and ~12 eV of the photo-absorption cross-sections cr(E) di- 
vided by the number of carbon atoms Nc, for the neutrals (dia- 
monds), cations (crosses) and dications (asterisks) of the PAHs 
considered, as a function of Nc- 



tral, cationic or dicationic are to account for the far-UV rise 
of the interstellar extinction curve they must also contribute to 
the short-wavelength side of the interstellar extinction bump at 
-2175 A (-5.7 eV). On the other hand, the onset of the far-UV 
rise produced by PAHs ought to depend on their charge state, 
therefore comparison with observed extinction curves could 
provide an observational handle to estimate the average charge 
state of interstellar PAHs. 
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The photo-stability of PAHs ++ deserves further investiga- 
tion: given the photo-absorption spectra we obtained, PAHs ++ 
will efficiently absorb UV photons, if available, and can be ex- 
pected to convert at least a fraction of the resulting energy to in- 
ternal excitation, possibly reaching their dissociation threshold. 
Should this be the case, this might turn out to be a major de- 
struction channel for PAHs. However, PAHs ++ absorb slightly 
less than their parent neutrals and singly ionised species be- 
tween ~1 and ~ 12 eV, and cool faster by IR fluorescence since, 
as we saw before, their vibrational transitions are overall from 
two to five times more intense. Since in these molecules uni- 
molecular dissociation is expected to occur in competition with 
vibrational cooling, PAHs ++ should actually be much more sta- 
ble against photodissociation than PAHs and PAHs + ;/ the dis- 
sociation thresholds are nearly unchanged by ionisation. More 
laboratory studies on the relaxation paths of PAHs ++ following 
photon absorption are needed. 



5. Conclusions and future work 

We performed a systematic theoretical study on the ionisation 
energies of PAHs and the IR and visible-UV absorption prop- 
erties of their dications. The values we found for the second IEs 
confirm that PAHs could reach the doubly-ionised state in Hi 
regions. Combining the IR and visible-UV absorption proper- 
ties of PAHs ++ we predict an increased photostability for these 
species with respect to their neutral and singly-ionized parent 
molecules, unless they have lowest dissociation channels. 

The recent paper bv lWittet alJ f2006) has recently renewed 
the interest in the dications of PAHs, proposing them as plausi- 
ble candidates of the Extended Red Emission. The fundamental 
data we present here for a large sample of different PAHs will 
be useful to model in detail the photophysics of these molecules 
in interstellar conditions. For instance, we can estimate quanti- 
tatively the complete IR emission spectrum of each molecule in 
the sample as a function of its charge state, for direct compari- 
son with data from the Infrared Space Observatory and Spitzer 
Space Telescope missions and, in th e near future, from th e 
Herschel Space Observatory mission (Mul as et alJl2006albld) . 
In addition, combined with measurements o f the electron 
recom bination rates of singly-ionized PAHs jNovotnv et alJ 
120051) . and their extension to doubly-ionized species, which 
still have to be studied, the present work opens the way 
for a detailed modelling of the che mistry and photophysics 
of PAHs in int erstellar environments (Ba kes & Tielensl ll998: 
iLePaee et all200ll) where they may play a prominent role. 
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Appendix A: Computed adiabatic and vertical 
single and double lEs 

The following Table. |A~T1 lists the single and double adiabatic 
and vertical IEs of each molecule considered in this work. As 
explained in the body of the paper, the single and double AIEs 
were computed, respectively, via total energy differences be- 
tween the geometry optimised cations and dications and the to- 
tal energy of the geometry optimised neutral. The correspond- 
ing VIEs were evaluated as the differences between the single 
point energy of the cation and dication at the optimised neu- 
tral geometry and the total energy of the optimised neutral. 
The adiabatic second ionisation energy Al is simply given by 
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Ka~Ka w hile its vertical value is obtained through the differ- 
ence between the single point energy of the dication at the op- 
timised cation geometry and the total energy of the optimised 
cation. We also evaluated the energy difference between the 
direct vertical double ionisation process (PAH— >PAH ++ ) and 
the intermediate one in which the mono-cation is allowed to 
relax to its ground-state geometry before second ionisation 
(PAH^PAH + ^PAH ++ ). The last two columns of Table |A~T1 
report the sum of the first and second vertical ionisation ener- 
gies /y + A/v and the difference between the double vertical 
ionisation energy 7+ + and + A/y 

For all singlet dications under study it is found that the in- 
termediate formation of the monocation is less energy demand- 
ing than the direct double ionisation process (the above differ- 
ence being positive); the opposite sign is instead ob tained for 
almost all triplet dications. As pointed out by ISchroder et alJ 
J200ll) . this comes from the triplet dication geometries being 
more similar to the neutral ground-state geometries than to the 
ones of the monocation. 



Table A. 1 : Adiabatic and vertical ionisation energies (in eV) of each PAH in our 
sample evaluated at the B3LYP/6-31+G* level of theory. Whenever available we 
list laboratory data for com parison. T he experimental single IEs come from the 
NIST Chemistry WebBook jLiasl2005l) while second and double IEs are, re spec- 
tively, from charge stripping and photon impa ct measurements oflSchroder et alj 
(2001) for corannulene and coronene and from lTobita e t al. 1 1994) in all the other 
cases. According to our calculations the dication electronic ground-states are 
marked in boldface. 



PAH 




First IE (/+) 


Dication 


Second IE (A/) 


Double IE (/ ++ ) 






molecule 


Adiab. 


Vert. 


Exp. 


state 


Adiab. 


Vert. 


Exp. 


Adiab. 


Vert. 


Exp. 


I* + A/ v 


/++ - (/+ + A/ v ) 


Azulene 
(C 10 H 8 ) 


7.15 


7.27 


7.42±0.02 


singlet 

triplet 


12.77 

13.39 


12.87 

13.53 




19.92 

20.54 


20.34 

20.63 




20.14 

20.80 


+ 0.20 

-0.17 


Naphthalene 
(C 10 H 8 ) 


7.80 


7.89 


8.144±0.001 


singlet 

triplet 


13.19 

13.65 


13.28 

13.77 


14.6+0.5 


20.99 

21.45 


21.35 

21.57 


21.5+0.2 


21.17 

21.66 


+ 0.18 

-0.09 


Acenaphtylene 
(C 12 H 8 ) 


7.66 


7.86 


8.12±0.10 


singlet 
triplet 


12.81 
12.78 


13.06 
12.88 




20.48 
20.45 


21.32 
20.66 


21.4+0.2 


20.92 
20.74 


+ 0.39 
-0.08 


Biphenylene 
(C 12 H 8 ) 


7.22 


7.37 


7.58±0.03 


singlet 

triplet 


12.09 

13.23 


12.26 

13.31 




19.31 

20.45 


19.94 

20.61 


19.7+0.2 


19.63 

20.68 


+ 0.31 

-0.08 


Acenaphthene 
(Ci2H 10 ) 


7.37 


7.46 


7.75±0.05 


singlet 

triplet 


12.42 

13.22 


12.51 

13.35 




19.78 

20.59 


20.17 

20.72 


21.7+0.2 


19.97 

20.81 


+ 0.20 

-0.09 


Fluorene 

(C13H10) 


7.56 


7.69 


7.91±0.02 


singlet 

triplet 


12.39 

12.86 


12.52 

12.99 




19.96 

20.43 


20.47 

20.59 


21.0+0.2 


20.22 

20.68 


+ 0.26 

-0.09 


Anthracene 
(C14H10) 


7.02 


7.09 


7.439±0.006 


singlet 

triplet 


11.68 

12.68 


11.74 

12.77 




18.70 

19.70 


18.95 

19.80 




18.83 

19.86 


+ 0.12 

-0.06 


Phenanthrene 

(C14H10) 


7.53 


7.63 


7 sol +n nni 


singlet 

triplet 


12.11 

12.19 


12.23 

12.32 




19.63 

19.72 


20.08 

19.86 




19.86 

19.96 


+ 0.22 

-0.09 


Pyrene 
(C16H10) 


7.07 


7.14 


7.426±0.001 


singlet 

triplet 


11.57 

12.28 


11.64 

12.39 


14.4±0.5 


18.63 

19.35 


18.92 

19.44 


19.8+0.2 


18.78 

19.53 


+ 0.14 

-0.09 


Tetracene 
(Ci 8 H 12 ) 


6.49 


6.55 


6.97±0.05 


singlet 

triplet 


10.66 

11.47 


10.70 

11.55 


13.1+0.5 


17.15 

17.96 


17.34 

18.11 


18.6+0.2 


17.25 

18.10 


+ 0.09 

+ 0.01 


Chrysene 
(Ci 8 Hi 2 ) 


7.17 


7.25 


7.60±0.01 


singlet 

triplet 


11.28 

11.67 


11.37 

11.75 




18.45 

18.84 


18.78 

18.94 




18.63 

19.01 


+ 0.16 

-0.06 


B enz [a] anthracene 
(Ci 8 Hi 2 ) 


6.98 


7.05 


7.45± 0.05 


singlet 

triplet 


11.16 

11.35 


11.24 

11.45 




18.14 

18.33 


18.42 

18.48 




18.29 

18.50 


+ 0.13 

-0.02 


Triphenylene 
(Ci 8 Hi 2 ) 


7.52 


7.86 


7.87±0.02 


singlet 
triplet 


11.72 
11.57 


11.81 
11.65 




19.24 
19.09 


19.59 
19.20 




19.67 
19.51 


-0.08 
-0.31 
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PAH 




First IE (/ + ) 


Dication 


Second IE (A/) 


Double IE (1 


++ ) 






molecule 


Adiab. 


vert. 


Exp. 


state 


Adiab. 


vert. 


Exp. 


Adiab. 


\Tart 

vert. 


Exp. 


/+ + A/ v 


/J + -(/+ + A/ v ) 


Corannulene 

(C20H10) 


7.62 


7.73 


7.83±0.02 


singlet 
triplet 


12.02 
11.66 


12.15 
11.76 


13.8±0.3 


19.63 
19.28 


20.14 
19.50 


20.1 ±0.2 


19.88 
19.49 


+ 0.26 
+ 0.01 


Benzo[a]pyrene 
(C20H12) 


6.71 


6.78 


7.12±0.01 


singlet 

triplet 


10.75 

11.47 


10.82 

11.58 




17.46 

18.18 


17.73 

18.30 


— 


17.60 

18.36 


+0.13 

-0.06 


Benzo[e]pyrene 

(C20H12) 


7.05 


7.12 


7.43±0.04 


singlet 

triplet 


11.20 

11.34 


11.26 

11.42 




18.25 

18.38 


18.51 

18.52 




18.38 

18.54 


+ 0.13 

-0.02 


Perylene 

(C20H12) 


6.57 


6.64 


6.960±0.001 


singlet 

triplet 


10.62 

12.03 


10.68 

12.11 




17.19 

18.60 


17.46 

18.73 


18.0±0.2 


17.33 

18.75 


+ 0.14 

-0.02 


Anthanthrene 


6.46 


6.52 


6.92±0.04 


singlet 


10.43 


10.48 




16.90 


17.10 




17.00 


+ 0.10 


(C 2 2H 12 ) 


triplet 


11.38 


11.54 




18.01 


18.09 


— 


18.06 


+0.03 


Benzo[g,h,i]perylene 

(C22H12) 


6.79 


6.86 


7.17±0.02 


singlet 

triplet 


10.74 

11.19 


10.80 

11.28 




17.53 

17.98 


17.78 

18.09 


— 


17.65 

18.14 


+ 0.13 

-0.05 


Pentacene 


6.12 


6.16 


6.63±0.05 


singlet 


9.91 


9.94 


12.8±0.5 


16.03 


16.18 


17.4+0.2 


16.11 


+ 0.07 


(C22H14) 


triplet 


10.56 


10.63 


16.67 


16.79 


16.79 


0.00 


Coronene 
(C24H12) 


7.02 


7.08 


7.29±0.03 


singlet 
triplet 


10.87 
10.72 


10.92 
10.79 


12.8+0.3 


17.88 
17.74 


18.11 
17.80 


18.5+0.2 


18.00 
17.87 


+ 0.11 
-0.06 


Dibenzo [b ,def] cry sene 
(C24H14) 


6.39 


6.45 


6.82 


singlet 

triplet 


10.11 

11.00 


10.17 

11.06 




16.50 

17.39 


16.73 

17.51 


— 


16.62 

17.51 


+ 0.11 

0.00 


Dibenzo[cd,lm]perylene 

(C26H14) 


6.31 


6.38 


6.72±0.02 


singlet 

triplet 


9.88 

11.03 


9.94 

11.12 




16.20 

17.35 


16.45 

17.43 




16.32 

17.50 


+ 0.13 

-0.07 


Hexacene 
(C26H16) 


5.83 


5.87 


6.36±0.02 


singlet 

triplet 


9.34 

9.84 


9.37 

9.90 




15.18 

15.68 


15.30 

15.78 


— 


15.24 

15.78 


+0.06 

0.00 


Bisanthene 


5.85 


5.90 


6.30 


singlet 


9.47 


9.51 




15.33 


15.48 




15.40 


+ 0.08 


(C28H14) 


triplet 


10.69 


10.76 




16.54 


16.70 


— 


16.66 


+ 0.04 


Benzo[a]coronene 

(C28H14) 


6.75 


6.81 


7.08 


singlet 

triplet 


10.34 

10.46 


10.39 

10.52 




17.10 

17.21 


17.30 

17.29 


— 


17.20 

17.33 


+ 0.10 

-0.04 


Dibenzo [be ,kl] coronene 
(C30H14) 


6.01 


6.05 


6.42±0.02 


singlet 

triplet 


9.54 

10.19 


9.57 

10.27 




15.55 

16.20 


15.70 

16.36 


— 


15.62 

16.32 


+ 0.08 

+ 0.04 


Dibenzo [be ,ef] coronene 
(C30H14) 


6.12 


6.17 


6.50 


singlet 

triplet 


9.66 

10.67 


9.69 

10.78 




15.78 

16.79 


15.94 

16.87 




15.86 

16.95 


+ 0.08 

-0.08 


Terrylene 

(C3qH 16 ) 


5.98 


6.05 


6.42+0.02 


singlet 

triplet 


9.33 

10.40 


9.38 

10.45 




15.31 

16.38 


15.55 

16.49 




15.43 

16.50 


+ 0.12 

-0.01 
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PAH 




First IE (/ + ) 


Dication 


Second IE (A/) 


Double IE (/ ++ ) 






molecule 


Adiab. 


vert. 


Exp. 


state 


Adiab. 


vert. 


Exp. 


Adiab. 


Vert. Exp. 


I* + A/ v 


/+ + -(/++A/ v ) 


Ovalene 


6.36 


6.41 


6.71 


singlet 


9.82 


9.86 




16.18 


16.34 


16.26 


+ 0.08 


(C32H14) 


triplet 


10.30 


10.35 




16.66 


16.72 ~~ 


16.76 


-0.04 


Tetrabenzo[bc,ef,kl,no]coronene 

(C36H 16 ) 


5.41 


5.44 




singlet 

triplet 


8.75 

9.66 


8.77 

9.72 




14.16 

15.07 


14.24 

15.17 ~~ 


14.21 

15.15 


+ 0.03 

+ 0.02 


Circumbiphenyl 

(C38H16) 


6.46 


6.52 


6.81±0.02 


singlet 

triplet 


9.58 

9.70 


9.64 

9.75 




16.04 

16.16 


16.26 

16.21 ~~ 


16.15 

16.27 


+ 0.11 

-0.06 


Circumanthracene 

(C 4 oH 16 ) 


5.90 


5.94 




singlet 

triplet 


9.10 

9.76 


9.13 

9.82 




15.00 

15.66 


15.12 

15.95 ~~ 


15.06 

15.75 


+ 0.06 

+ 0.20 


Quaterrylene 

(C40H20) 


5.62 


5.68 


6.11±0.02 


singlet 

triplet 


8.50 

9.36 


8.55 

9.40 




14.12 

14.97 


14.35 

15.08 ~~ 


14.23 

15.08 


+ 0.12 

0.00 


Circumpyrene 

(C42H 16 ) 


6.06 


6.10 




singlet 

triplet 


9.13 

9.68 


9.16 

9.73 




15.19 

15.75 


15.35 

15.80 


15.27 

15.84 


+ 0.08 

-0.04 


Hexabenzocoronene 


6.55 


6.79 




singlet 


9.58 


9.63 




16.14 


16.33 


16.42 


-0.09 


(C42H18) 




triplet 


9.52 


9.57 




16.07 


16.15 


16.36 


-0.21 


Dicoronylene 

(C48H20) 


6.07 


6.13 




singlet 

triplet 


8.84 

9.37 


8.89 

9.40 




14.91 

15.44 


15.13 

15.50 ~~ 


15.02 

15.53 


+ 0.11 

-0.03 


Pentarylene 

(C50H24) 


5.37 


j.hj 




singlet 

triplet 


7.92 

8.64 


7.97 

8.68 




13.29 

14.00 


13.51 

14.10 


13.39 

14.11 


+ 0.11 

-0.01 


Circumcoronene 


6.14 


6.35 




singlet 


8.90 


9.04 




15.05 


15.19 


15.28 


-0.09 


(C54H18) 




triplet 


8.81 


8.85 




14.95 


15.00 


15.19 


-0.19 


Circumovalene 
(C66H20) 


5.71 


5.74 




singlet 

triplet 


8.32 

8.67 


8.34 

8.70 




14.02 

14.37 


14.12 

14.42 


14.07 
14.44 


+0.05 
-0.02 



